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Abstract 



We apply the gaussian variational method (GVM) to study the equilibrium 
statistical mechanics of the two related systems: (i) classical elastic manifolds, 
such as flux lattices, in presence of columnar disorder correlated along the r 
direction (ii) interacting quantum particles in a static random potential. We 
find localization by disorder, the localized phase being described by a replica 
symmetry broken solution confined to the mode cj = 0. For classical systems 
we compute the correlation function of relative displacements. In o? = 2 + 1, 
in the absence of dislocations, the GVM allows to describes the Bose glass 
phase. Along the columns the displacements saturate at a length /_|_ indicating 
flux-line localization. Perpendicularly to the columns long range order is 
destroyed. We find divergent tilt modulus C44 = 00 and a x ~ r^/^ scaling. 
Quantum systems are studied using the analytic continuation from imaginary 
to real time r — > it. We compute the conductivity and find that it behaves 
at small frequency as (7{uj) ~ uj"^ in all dimensions {d < 4) for which disorder 
is relevant. We compute the quantum localization length ^. In d = 1, where 
the model also describes interacting fermions in a static random potential, we 
find a delocalization transition and obtain analytically both the low and high 
frequency behavior of the conductivity for any value of the interaction. We 
show that the marginality condition appears as the condition to obtain the 
correct physical behavior. Agreement with renormalization group results is 
found whenever it can be compared. 



I. INTRODUCTION 



Two related longstanding problems of statistical mechanics have recently received re- 
newed attention. On the one hand the problem of the pinning of the flux lattice in high-Tc su- 
perconductors has motivated further studies of the effect of disorder on an elastic manifold§. 
Pinning by correlated disorder, e.g. in the form of columnar defects introduced in the sample 
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by heavy-ion irradiationi'i, is particularly interesting for technological applications^^ since 
it leads to a large increase in critical current. On the other hand, the problem of inter- 
acting quantum particles in a random potential, as the one of quantum fluctuations of an 
elastic manifold in a random potential, have been recently studied by a number of authors. 
The latter problems have many experimental realizations, such as Charge Density Waves 
Wigner crystalJli, quantum vortex creeD in superconductorsB, superconducting- 
insulating transition in disordered superconductorsliy. There is a deep analogy between these 
two classes of problems which rests on the well known identification of the bosons worldlines 
in the imaginary-time path integral of quantum mechanics with actual flux lines0. This 
relation has been reanalyzed&B'l in details and has led for instance, to the prediction of 
a new physical state of vortex system^ with vanishing linear resistivity, where 3d flux lines 
are localized by columnar disorder, very much like 2d quantum bosons become localized 
in presence of disorder and form the so-called Bose glass statdll. Although there IS some 
evidence for a Bose glass, or localized phase, few analytical methods are available to study 
this problem. One has to resort to series expansions, numerical simulationsEHnl. Exper- 
imental evidence is also controversial. In the case of d = 1 (one spatial dimension, one 
"time" dimension) where interacting bosons and fermions can be related, the problem can 
be studied using powerful bosonization and Renormalization Group (RG) techniques. These 
methods showed a transition away from the superfiuid phase when the repulsion is increased 
and allow to compute the crossover towards the localized phase. However they do not allow 
to describe the low energy properties of the resulting (presumably) localized phase, since 
RG trajectories runaway to strong coupling. Only in the case of non-interacting fermions, 
or equivalently bosons interacting with hard-core repulsion, is the "Bose glass" phase rig- 
orously known to exist. The conductivity was calculated in that case by Berezinskii@ and 
found to behave as <j{uj) ~ cu^ ln(c<j)^. 

An interesting method to study problems with an interplay of elasticity and disorder 
is the Gaussian Variational method (GVM) introduced by Mezard and Parisiii to study 
the problem of elastic manifolds, such as interfaces, in a random uncorrelated potential. 
In this method one finds the best quadratic Hamiltonian approximation to the original 
complicated Hamiltonian. This quadratic approximation is performed in replica space and 
the optimization naturally leads to Replica Symmetry Breaking (RSB) whenever the energy 
landscape is complicated and contains metastable states. Within the gaussian ansatz, RSB 
appears as the natural way to encode the distribution of the metastable states induced 
by disorder. The possible optimal Hamiltonian is parametrized by a function a{u), which 
leads to a large number of optimization parameters. Thus one can hope that it produces a 
reasonable approximation of the real problem. Indeed, in the interface problem, the method 
leads to a reasonable, Flory-like, approximation for the roughness exponent; It also leads 
to the correct physics of sample to sample fiuctuationsil in response functions, due to rare 
events. This was recently demonstrated in d = 1 -|- 1 using the mapping onto the Burgers 
equation^. In recent papersilS, we have applied the GVM to the problem of an elastic 
lattice in a random uncorrelated potential. It was shown&^i that the periodicity of the 
lattice leads to logarithmic disorder-induced displacement at large scales in d > 2. 
Comparison with two RG calculationsc3'ElS suggests that the GVM captures the correct 
physics, and leads to reasonable quantitative approximations. Although the GVM captures 
important non perturbative features such as the existence of many local equilibrium states, 
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it may miss some details of the nonlinearities and instanton-type configurations (kinks) 
between metastable states. 

In the present paper we apply the Gaussian Variational Method to classical problems with 
correlated disorder and to quantum problems with disorder. A summary of the results of this 
work has appeared in 0. There are some additional subtleties in applying the method, due 
to translational invariance in one direction ( the "imaginary time" direction) and extra care 
should be paid in treating boundary conditions and taking the large "imaginary time" limit 
/3 — > cxD properly. We then find that the GVM naturally leads to localization by disorder 
and we compute the properties in the localized phase. For the classical system we obtain 
the tilt modulus and find that 044(5^) — > 00 as — > in the localized phase, i.e a signature 
of the transverse Meissner effect predicted in Ref. |. For quantum systems we investigate 
in detail dynamical correlation functions such as the conductivity. We study in detail the 
case of of (i = 1 which describes either the classical problem of flux lines in a plane with 
correlated (columnar) disorder or interacting fermions (or bosons) in a d = 1 static random 
potential. We also give the main features and the essential physics of the solution for general 
d. In (i = 2 we study simplified models relevant to experimental systems such as the three 
dimensional vortex lattice with columnar disorder, disordered quantum bosons in c? = 2, the 
pinned Wigner crystal and Charge Density waves. 

In a celebrated workS, Fukuyama and Lee (FL) analyzed the quantum problem of a 
one-dimensional CDW with disorder [d = 1) using a diagrammatic self-consistent method. 
Physically, their method, and its further applicationsS to interacting fermions m d = 1, 
amounts to choose an unknown classical solution {h = 0) of the equations of motion in the 
random potential and to treat the quantum fluctuations as an harmonic oscillator around 
this solution. Since it corresponds to an expansion around an unknown classical background, 
the FL method contains an ad-hoc phenomenological parameter. Although the static cor- 
relations are nearly independent of this ad-hoc parameter, completely different physical be- 
haviors can be obtained for frequency-dependent quantities depending on its precise value. 
To obtain the presumed correct physics the parameter had therefore to be adjusted by hand. 
Remarkably, the physics of the GVM solution parallels the FL theory and allows to estab- 
lish it on a firmer footing. As in FL, the physics of the GVM solution involves an elastic 
expansion around the classical solution, but this comes out naturally as the solution of the 
variational equations. In addition the GVM allows for a detailed description of the classical 
solution. In particular, the phenomenological parameter adjusted by hand in FL exactly 
corresponds to the breakpoint parameter Uc in the RSB solution ! It is thus determined by 
the theory itself. The GVM therefore provides a more general theoretical framework which 
can easily be extended to higher dimensions and more complicated systems. 

In the models considered here, weak pinning is assumed (disorder is weak and gaussian). 
We have not attempted to treat the strong pinning case, i.e a Poisson distribution of strong 
pins, although this extension of the method is possible in principle. For flux line lattices, 
disorder can be considered as gaussian and weak when > D where D is the distance be- 
tween columnar defects and the Larkin-Ovchinnikov pinning lengthiHi. Since columnar 
defects produced by irradiation are quite strong pinning centers this corresponds to high 
fields or temperature. Furthermore our theory corresponds to a large pinning length ^ ^ a, 
i.e the elastic limit is assumed, as well as the absence of dislocations. 

The plan of the paper is as follows: in Section y we introduce the two classes of models 
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studied in this paper, quantum problems and classical problems with correlated disorder. In 



Section |ITl| we apply the variational method to these models and derive the corresponding 
saddle point equations. We discuss general features of the solution and show that, within 
the GVM, only the zero Matsubara frequency mode can exhibit replica symmetry breaking. 
In Section |IV| we give the solution for a classical model in ci = 1 and d = 2 describing respec- 
tively lines and vortices in a correlated random potential, as well a the solution for general d. 
In Section |V| we discuss the properties of the solution specific to quantum problems, such as 
the analytical continuation to real time and we perform the calculation of the conductivity. 
The GVM is applied to the problem of interacting quantum particles in d=l+l. Conclusions 
can be found in section [V^. 

II. MODELS 

This section introduces the various classical and quantum models which we study using 
the variational method. 



A. general model 

A general model for an elastic manifold of internal dimension d+1 in a correlated random 
potential is: 

H = J d'^xdT^[{V^(j)f + {dr<pf] + J d'^xdTW{{(f){x,T),x) (l) 

(/) can in general be a A^-component field (to describe a manifold embedded in N dimensions) . 
The partition function is Z = J Z)0(x, r) exp(— ■^). The disorder potential is Gaussian, 
correlated in the direction r and the correlator is: 



W{^, x)W{{^', x') = - ^')6\x - x') (2) 

This model describes a large number of classical models and we give important physical 
realizations in the next sections. It also describes the equilibrium statistical mechanics of 
quantum problems with disorder. The Hamiltonian (P can also be viewed as the action of 
a quantum problem and T corresponds to h. t = it is the imaginary time of the quantum 
problem in the Matsubara representation. (3 is the inverse temperature of the quantum 
problem, and by imposing periodic boundary conditions (j){x, t = L) = (f){x, r = 0) one 
describes bosonic systems, with L = (3h. Disorder is then time independent which is the 
natural situation for a quantum problem, d = corresponds to a single quantum particle 
in a random potential. Applications to quantum problems are considered in more details in 
Section ^ Note that in this paper [3 and T are two unrelated parameters. 

In order to average over disorder one uses the replica trick by introducing k replicas of 
the system 0^, a = 1, ..k and taking the limit — at the end. The replicated Hamiltonian 
averaged over disorder reads: 

^ = J d'xdr^ E[(V.0a)' + (9.0.)^] + d'xdrdr' ^ r) - r')) (3) 

1 a,b 
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B. vortex lattice: all harmonics 



An important physical realization of the system is the Abrikosov vortex lattice in 
the presence of columnar disorder. The columns are oriented parallel to the r direction, the 
magnetic field being aligned with the columns. The d = 1 version of this problem describes 
flux lines in a plane pinned by columns of defects while the d = 2 version describes the flux 
lattice with columns aligned with the direction of the flux lines. 

We denote by Ri the equilibrium position of the vortex lines labeled by an integer i, 
which defines a perfect lattice of spacing a, and by u{Ri, r) their displacements transverse 
to the field. As will be seen later, the disorder defines a pinning length ^ (in the transverse 
direction) such that the relative displacement of vortices separated by ^ (for equal r) is of 
order a. For weak disorder a/^ -C 1 it is legitimate to assume that u{Ri, r) is slowly varying 
on the scale of the lattice and to use a continuum elastic energy, in terms of the continuous 
variable m(x, r). Impurity disorder is modeled by a r independent gaussian random potential 
V{x) with correlations: V {x)V [x') = W6{x — x'). The total energy is: 

H=^J d^xdr c[{druf + {V ^uf] + j dxdTV{x)p{x, r) (4) 

where the density is p{x,t) = J^i^i^ — Ri — u{Ri,T)). For simplicity we use an isotropic 
elastic hamiltonian. Instead of using the displacements u, it is more useful to introduce a 
smooth "labelling" fieldS^ ri{x, r) = x — u{ri{x, t), r). The density can be rewritten as 

pix, t) = podet[d^rjp] ^ e^^''^^'^) ^ Po(l - d^uM^, r), r) + ^ e^^>^(x)) (5) 

where the usual translational order parameter defined in terms of 

the reciprocal lattice vectors K and po is the average vortex density. 
Using the replica trick on (11) the disorder term gives 



W_ 
2T 



J dxdrdr' p°'{x,t)p^{x,t') (6) 

a,b 



The above decomposition for the density leads to: 

Hes=\l d'^xdr ^ ciidru'^f + (V.n")2] (7) 

[^d^U%X,T)dyix,T')+j:^'^ 
a,b '^^ K^O 



J d''xdTdT'J2[^d^u''ix,T)d^u\x,T') + J2 ^COsiKiu''ix,T)-u\x,T' 



The range of validity of (0) is as discussed in Ref. ^,^: (i) elastic limit a/^ <^ 1 (ii) neglect 
of dislocations, (which is fully justified in = 1). We redefine p^W ^ W in the following. 

The slowly varying part of the disorder, i.e the g ~ Fourier components of V{q), couples 
to the long wavelength part of the density fiuctuations and gives rise in (^ to the quadratic 
part. The higher Fourier components of the random potential, specifically those of Fourier 
components near q ^ K, lead to the cosine terms. The g ~ part of the disorder can easily 
be treated, either in the replica form of (|^, or before averaging by simply redefining: 
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Ua{x,T) Uaix,T) + fa{x,T) fa{q,Uj) = S^fi " (8) 

where V^^o is a truncation of the random potential V{q) keeping only the Fourier components 
close to zero. Since this redefinition of the u field is r- independent, it does not affect 
the cosine terms in The g ~ part of the disorder decouples completely from the 
Hamiltonian (0) and we thus drop it in the remainder of this paper. It leads to a simple 
additional contribution to the correlation functions which can be computed using (||). 

The model which describes vortex lattices in correlated disorder is thus a particular 
case of the general model (|l],|^. 



C. single cosine model 

A further simplification can be obtained by keeping only the lowest harmonic Kq in (^. 
This leads to another model, which for simplicity we define for a scalar field (p. This model 
is studied in detail in this paper: 



^'"^ = J d'xdr^ E(V.0a)' + {d^' - J d'^xdrdr' ^ cos(2(0,(x, r) - r'))) (9) 

a,b 



It applies directly to CDW and interacting fermions in li = 1 as discussed in Section 0. It 
is also relevant for fiux lines and correspondence with model (|^ can be made by defining 
u = atp/n. As we will see, in the limit of weak disorder a/^ <^ 1, the contribution of higher 
harmonics becomes irrelevant at large scales, for separations transverse to the magnetic field. 
As for point disorder, an intermediate "random manifold" regime exists at intermediate scale 
X < ^ (see Ref. ^9[). For separations along r, higher harmonics can be included, as in (^, 
but do not change the main physics. 



D. domain of application of the models 

In this paper we consider Gaussian disorder, i.e a large number of very weak pins, rather 
than the Poissonian disorder which corresponds to few strong pins. In addition, we are in 
the elastic limit of small relative displacements between nearest neighbors. Finally in ci = 2, 
dislocations have been neglected. 

Because of the elastic forces (interactions) the modelization by Gaussian disorder is often 
appropriate. Its domain of validity can be estimated at low T as follows. Let us consider 
defects of strength Uq separated by a distance i.e a density n ~ 1/D'^. The disorder 
potential is V{x) = UoJ2i^{x — Xj) where Xj are the positions of the defects. Using (^, 
and keeping the lowest harmonic, the coupling to density is ~ pQj2iCos{Ko.{u — Xi)). If 
density of defects is very small, or Uq very strong, the phase of the lattice adjusts such that 
Ko.u = Ko.Xi at each defect. In fact one has to balance the cost in elastic energy with the 
gain in potential energy. Assuming that Kq.u varies by 2tt over distances of order ^, the 
elastic (kinetic) energy is c ^"^"^ and the potential energy is, assuming effectively a Gaussian 

disorder, UoJn^'^. Optimizing over ^, one finds ^ ~ (l/nUo'^y^^^''^^ The Gaussian-like 



6 



regime fails when the distance ^ becomes smaller than the distance D between defects. One 
then enters the Poissonian regime where ^ then saturates at ^ = D. From the expression 
$^/D = (n^(2-d)/rf^^^2-ji/(4-d) Qj^g finds the domain of vahdity of Gaussian disorder in any 
dimensions. In c? = 1 one has ^/D = (n/Uoy^^ and the Gaussian regime holds if the density 
of defects is not too small n ^ Uq. In d = 2 one has C,/D = I/Uq and the Gaussian regime 
holds at any density of defects provided disorder is small enough f/o <^ 1. In d = 3 one 
has ^/D = 1/ [n^/^Uo^] and the Gaussian regime holds for all density and disorder provided 
they are small. 

For models described by a single harmonic, such as charge density waves, the above 
argument always holds. For flux lattices, disorder can vary at scales much smaller than a 
and higher harmonics may be importantlH. There is an additional length scale, the Larkin 
Ovchinnikov length i?c, for which relative displacements are of order ^o- At weak disorder 
Rc and ^ (also called Rg) are very different^. A sufficient criterion for Gaussian disorder to 
hold is then that Rc^ D. 



III. VARIATIONAL METHOD : GENERAL PROPERTIES 

A. saddle point equations 

We study the statistical mechanics described by the partition function: 

D0(x,r)exp(-^) (10) 

for a system of finite size j3 in the r direction. H is the hamiltonian (|l[). In view of further 
applications to quantum problems we impose periodic boundary conditions in the r direction. 
For infinitely thick vortex systems of length L = (3 ^ oo boundary conditions should not 
matter. However, imposing definite boundary conditions and working first on a finite size 
sample is crucial for obtaining well-defined equations in the case of correlated disorder. 
Introducing replicas and averaging over disorder the model becomes (^). We now use 



the variational method as in Ref. ^ This is done by choosing a Gaussian variational 
Hamiltonian: 

1 f d'^q 

^0 = - y j^J2'^abil,^n)Ml,^n)M-1,-^n) (H) 

where the discrete summation is over Un = 2im/ (3 analogous to Matsubara frequencies. Note 
that ^ Eg ^ / (0r and ^ Ea;„ ^ / |^ in the large volume (fi, resp. (3) limit. 

By the same arguments as in Ref. ^ this method is exact for a choice of correlator of the 
form C{(j) — (f)') = NV{{(f) — (p'Y/N) in the limit oo. This corresponds to a manifold 

in infinite dimensions, either a classical manifold awith correlated disorder or a quantum 
manifold with static disorder. We obtain the variational free energy Fvar = Fo+ {H — Hq)ho 
as 

^var = I E E + ?')G(g, U^a ' [In(TG)], 

+ E / f f drdr' V{Babix = 0, r - r')) (12) 
21 ^ J Jo Jo 
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where 

Bai>{x,r) = ([0,(x,r) -0,(O,O)]2) (13) 

= -B I 7 7d ^{Gaaiq, UJn) + Gbb{,q, UJn) - 2 COs(gX + UnT)Gab{q, UJn)) 

In general the function V{B) which appears in the variational energy is defined as: 

{cm = vm) (14) 

where (..) denotes an average over any (one-component) gaussian random variable 0. For 
the case of model (g) V{B) = -W exp{-2B) and V{B) = -VT Ei^ exp -^2^/2 for model 
(^). Defining the self-energy = cq^5ab — <^ab-, and G^^{q) = Y.bGab{q) we obtain by 
minimization of the variational free energy Fvar the saddle point equations: 



G-' = c{q' + ujD + lf dr{l - cosM)[V'{B,M r)) + E ^'(^^.(O, r))] 

T ^0 b^a 

CTa^b = — / rfrcos(u;„r)\/'(5„b(0,r)) (15) 
1 Jo 

Taking the limit of the number of replicas k ^ and considering the most general solution, 
i.e. a replica symmetry broken solution, we denote G(g,to'„) = Gaa{q,^^n) and parameter- 
ize Gab{,q,^n) by G{q,Un,u) where < m < 1, and similarly for Bab{x,T) and B{x,t,u). 
Equ. ([1^) and the algebraic rules for hierarchical matricesH^ give: 

G:' = c{q' + u;l) + ^ f dr{l-cos{uJnr))[V\B{r))- [' duV'{B{r,u))] 
I Jo Jo 

2 

a{ujn,u) = — dTcos{uJnT)V'{B{T,u)) (16) 



T 



with 



2T f d'^q ~ 
B{r, u) = ^n) - G{q, Un, u) cos(w„r)] 



2T r/S d'^n ~ 

B{r) = y EX ^dGi^^^n)il - cos{u;^t)) (17) 

where we denote B{t) = B{x = 0,r) and B{t,u) = B{x = 0,t,u). A replica symmetric 
solution would correspond quantities constant in u in ([l6|) , and is analyzed in Appendix 
Note that the connected part G~^{q,uj), is renormalized by correlated disorder, while it is 
not the case for point disorderly. This is a crucial difference which leads to localization and 
stiffening of the tilt modulus. 



B. general features of the solution 

A general and important property of quantum problems with disorder is that off-diagonal 
quantities like Ba^bix, r) are in fact always independent of r. The general argument^a is 
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that in each reahzation of the random potential W, the disorder does not depend on r. 
Therefore before averaging over disorder: 



Gab,w = {Mx, r)0b(O, 0)) = r)) (0b(O, 0)) = {M^, 0)) (0,(0, 0)) (18) 

It is important to note that such a property crucially depends on the assumption that the 
hamiltonian is r-independent and on the fact that equilibrium has being attained. This is 
the case considered here. Perturbations which explicitly depend on imaginary or real time 



destroy this property and may lead to more complex behavior, as discussed in Section ^1. 

It is easy to check on the saddle point equations (|r^) that they indeed admit a solu- 
tion such that off-diagonal quantities are r-independent. For this solution the saddle point 
equations simplify greatly: 

G:' = c{q' + ujl) + ^ f dr{l-cos{u:nr))[V\B{r))- [' duV'{B{u))] 
1 Jo Jo 

23 

a{uJr,,u) = ^V'{B{u))6n,o (19) 

with 

Bin) = f E / + f / ^M^^^n = 0) - G{g,u^ = 0,u)] 

Bir) = —Y.l ^,GM,^n){l-cosM) (20) 

where we denote B{t,u) = B{u). The simplifications are the following. Since B{t,u) is 
independent of r the self energy cr(c<j„,'u) vanishes for all modes excepted Un = 0. One also 
has G{q, Un, u) = for c<j„ 7^ 0, and thus G{q, uJn) = Gc{q, uJn) for Un 7^ 0. 

An obvious property of the above saddle point is thus that replica symmetry breaking is 
confined to the mode Un = 0. The equation for this mode: 

2/? 

a{un,u) = 6n,o(^{u) = —6n,oV'{B{u)) (21) 



is identical to the one for a model with pointlike disorder in d dimensions, studied 
The only (important) difference is that the potential is effectively multiplied by f3, i.e the 
total length along the r direction. This is why a single mode {un = 0) can lead to finite 
contributions to correlation functions in the limit /3 — ^ 00. In that limit the corresponding 
model with point disorder in d dimensions can, in some sense, be considered as effectively 
at zero temperature ~ T/p. 

Known resultsil'ilii can thus be transposed to study the zero mode Un = 0. For the 
potentials with power law correlators!! V{x) = 5'x^~'^/(2(l — 7)) there are two generic cases: 
(i) long range potentials, if 7(1 — d/2) < 1, for which one has full replica symmetry breaking 
and (ii) short range potentials, if 7(1 — (i/2) > 1, for which one has one step replica symmetry 
breaking (and a transition to a high temperature Replica Symmetric (RS) phase). In the 
case of the single cosine model (|^) there iscan^ a one step solution for d < 2 and a full RSB 
solution for d > 2. Potentials V{x) which are realistic to describe media with several scales, 
such as the vortex lattice (|^), will include several of the above regimesSSS. 
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The general way to study ( pi]) is as follows. There is a breakpoint Uc such as (t{u) is 



constant for u > Uc- For u > Uc the inversion rules applied to (20) give: 



2T r rl'^n 2T r rl'^n 1 

Bin > u,) ^ BM . + -/ (22) 

where by definition Si = [a"](Mc) and the function [o"](m) = ua{u) — Jq dva{v). One searches 
for a full RSB solution by applying to (^) the inversion formula for u < Uc- 

S M = S Me + -T / dv — ^ (23) 

Differentiating ( pT]) with respect to u, using (|23| ) and [o"]'(m) = ua'{u) one can then divide 
by (t'{u) 7^ since full RSB is assumed. The quantities B{u), [o]{u) and a{u) are then 
determined by elimination from the resulting system: 



AV"{B{u)) 



{27^Y{cq^ + [a]{u)f 



a{u) = ^V'iBiu)) (24) 

and the definition of [cr](M). These equations, which govern the mode a;„ = 0, are identical to 
the one obtained previously for pointlike disorder in d total dimension, up to a rescaling of u 
by i.e u is smaller by a factor of (3 in this problem, [a] and B do not scale but a is larger by 
a factor p. The full RSB solution for the power law potential is thus [(t](m) = T,i{Pu/ucY^^ 
with 9 = {2 + {d — 2)7)/(l + 7). Once the equation is solved for u < Uc, and provided that 
a full RSB solution exists, which is the case when ^ > 0, the unknown constants Si and 
B{uc) are determined by matching or equivalently by elimination in: 



-W"{B{u,)) 



d'^q 



{2t:Y {cq^ + Si)2 
^ , , 2T f d'^q . . 2T f d'^q 1 

For the single cosine model (|^) which we focus on here, the above system ( p^) simplifies 
into: 

^ = ^"^")/(^(cg^+Un))^ (2^) 



and the solution of (|2|) isBESl in ci > 2 [a]{u) = {iiuju^fl^ and = 8Tqc-'^/V(4 - d). 
The exponent 9 = d — 2 characterize the scaling of the fluctuations of the energy of the 
metastable states with distance, AF ~ x^. A pair of states characterized by u in Parisi's 
hierarchy are typically x apart in space and AF = T/u apart in energy. The large scale 
behavior is thus controlled by small u. 

In c? < 2, the single cosine model does not admit a full RSB solution. Instead one finds 
a one step RSB solution with a{u) = for < and a{u) = cr{uc) for u > Uc- One then 
has: 



10 



[a]{u)=0 u<Uc (27) 
[a]{u) = ^, = Pu,^V'{B) u>u, (28) 

where we have used ( pi]) and the fact that Ei = [o"](uc) = UcCr(uc) for such a solution. 
One has also B{u) = oo for u < Uc and we denote B{u > Uc) = B. For the single cosine 
model El = /3uc^e~^^. To determine Si, B and Uc in the one step case, one can use only 
equations (p8| , |23|) . One equation is still missing. Before we discuss the additional equation 
needed to determine Uc, let us further simplify the saddle point equations. 

Irrespective of whether the solution is one step or full RSB, one can always rewrite the 
variational equations as: 

G;\q, ^„) = c(g2 + ujD + Si(l - 5„,o) + IM (29) 

IM = I f\r{l - cos(^„r))(y'(5(r)) - V\B)) (30) 
1 Jo 



with 



^ 2r /• A v-^/ N 2T r d^q i 

^ = lfy (2^5Go(,.c.„) + yy P^^j:^^ (31) 

2T cf^q 

where B{u > Uc) = B. To obtain (^) from (p!6D we have splitted: 

/ duV'{B{u)) = (1 - u^)V'{B) + / duV'{B{u)) = V'{B) - —Si (32) 
Jo Jo 2p 

and used the equation (^Tj) and the definition of [o"](m). The form (^) is convenient because 
when P becomes large, B{t) converges at large r towards B and thus /(u;„) in (|30D goes to 
zero at small One has: 

B - B{r) = =^ V / 7^ cos(cu„r)— — ^- — - (33) 

and thus limT-^oo lim/3^cxD B{t) = B. 

One thus sees on (^) that disorder produces a "mass" term Si and thus localization of 
the elastic manifold. The mechanism by which it is generated is subtle, however. Such term 
cannot exist at ci;„ = as can be seen from (|16D . Indeed G'^^{q = 0, c<j = 0) = as it must 
from translational invariance after averaging over disorder. As soon as Un = 27m/ P ^ the 
first term in (pOf) becomes non zero equal to Si while the second one is a smooth function 
of ujn and therefore remains small near u = 0. Thus we find, and observe in a numerical 
solution of equations (^^T]), that in the limit /5 — > cxd, G~^{q,uj) develops a discontinuity 



at u; = 0. In fact the only way to obtain a finite localization length in the limit /5 — > cxo 
is to break replica symmetry. The formula (0) for B shows that Si plays the role of an 
effective mass which makes B finite in d < 2. There is no such mechanism in the replica 
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symmetric solution, which cannot describe properly the localized phase and is studied in 
detail in Appendix ^ 

Remarkably, Eqs. (p9|-pl|) now forms a closed system of equations, depending only on 
one additional quantity Si, and do not depend on how the replica symmetry is broken or on 
the detailed structure of the solution B{u) for u < Uc- We will discuss further this property 



in Section VI 



IV. ELASTIC MANIFOLDS IN CORRELATED POTENTIALS 



It is possible to solve analytically the saddle point equations (p9| - |3TD in the limit of low 
temperature T — > 0. Such solution, as we will show, contains the essential physics of the 
localized phase. It is important to retain the finite temperature only for few quantities that 
we will discuss explicitly. 

In the limit when T — > 0, one can expand the potential in the expression ( PDD of I{uJn)- 
One gets a self-consistent equation for /(cu^): 

- -4V"(0) / - e,^ + S.+U + 

It is important to note that the low-T expansion of the potential can be performed safely 
only on the form (^) of the saddle point equations. Si being fixed. It would be incorrect to 
expand directly (|19D since B{u) diverges at small u and the low-T expansion breaks down. 
We now study different cases. 



A. Solution for the cosine model in d = 2 

In (i = 2 the solution for the single cosine model (|^) is one step but marginally so. 
The value of Uc is determined unambiguously either from the full RSB solution in the limit 
d ^ 2^, or by minimizing the free energy Fvar in (HID. This is done similarly to the case of 
point like disorderii. For model (^ in the continuous limit: 

and thus I3uc = T/vrc for model (P). Note that there is a transition to a high temperature 
phase where disorder is irrelevant at T = = 7t(3c. For f3 not very small this temperature 
however is very high and in realistic systems other mechanisms, e.g dislocation induced 
melting, will cause a transition out of the glass phase before this one occurs. Thus we 
confine our study to the low temperature glass phase T ^ Tc. 
One then obtains from 



IM = 4ncS, I - ^ ^ ^ (36) 



which leads to, in the limit where the cutoff goes to infinity 
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I{ujn) = Si log(l + ) (37) 

^1 



The solution of (|37|) can be written as 

JK) = J:^f{V^\uJ^\/J^^) (38) 



Si defines a natural length scale ^ = yc/Si. ^ is called the pinning length and corresponds 
to the scale above which the long range order is destroyed due to disorder. 
For small Un the solution of (137|) is 



I{uJn) ~ y2cSi|u;„| (39) 

and for large a;„ 

J(u;„)~Silog(ccu2/Si) (40) 
This formula is valid up to the cutoff ujn ~ A. At large ui, I{uj) — > /(oo) = / ^^^^ ^^^^^^ . Si 



is determined from (p5|,|35|,|28D . At low T one finds Si = AW/n. 



The general expression for the correlation functions for a one step solution is: 

+ 7^ / , , ^ , 1 - cos gx 41 

(5ucJ (27r)^cg2(cg2 + Si)^ ^ ' 

At X = 0, -8(0, r) goes to a constant when r ^ oo. This describes the localization of a flux 
line due to the correlated disorder. One defines a localization length l± as: 



a 



2 



^i = ^Jl^(['^(^.^)-<^(o.o)F) (42) 



remembering that we have defined the displacements as m = a(j)/T:. One finds from (^if) : 



^ 7r2 7r2 /3 (27r)^c^2+cg2 + Si + /(^„) ^ ^ 



In d = 2 one finds: 



,2 ,2 a2 2r /s7 1 r , , 1 + ,,,, 

where = 2Ta^/(c7r^) / (g^ + cu^)"^ is the Lindemman length. This formula can be 
used to compute the upward shift in the melting temperature due to correlated disorder 
using a Lindemman criterion. 

On the other hand for fixed r, the Un = mode dominates at large x and B{t, x) grows 

as: 
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r) ~ log(x/0 = log(^) (45) 

The prefactor in (^3|) is universal and lead to a power law decay of the translational order 
paramete li. The additional contribution of the g ~ Fourier components of the random 
potential can be computed from (^ and gives a small non universal correction. 

One can define the cu-dependent (i.e dependent, in usual notations) tilt modulus C44{uj): 

T 

\im{dr(f)dr(l)) = \irauj'^{(f){q,uj)(f){—q, —to)) (46) 



Using equation (|29|) one finds: 

1 a;2 



044(00) cuj'^ + Si + I{uj) 



(47) 



a formula which holds in any dimension. Thus we find that C44 — C44 (dJ 0) diverges in the 
Bose glass phase, a signature of the transverse Meissner effecta. 

We find from the GVM three important properties of the Bose Glass phase: the finite 
localization length l±, the decay of the translational order in the direction transverse to the 
columnar defects, and the divergence of the tilt modulus. 



B. Solution for the cosine model in d = 1 

In = 1 thermal fluctuations are much stronger and in the absence of disorder correlation 
functions grow logarithmically at large distances: 

Bix,T)r^ ^Jn[ix' + r')/a'] (48) 

In the presence of disorder one expects a low temperature localized phase were disorder is 
relevant and a replica symmetric high temperature phase were thermal fluctuations domi- 
nate. The transition at Tbg = Qttc/Kq (at infinitesimal disorder) can be seen on the RS 
solution (see appendix ^ and is found to coincide with the resultsii of the RG. Eq. (|34D 
leads to: 

/K) = -^^i^l'^' - (^n + Si + HoOn))-'/'] (49) 

To solve this equation one needs to know Si, and thus determine Uc- Most of the physical 
properties, such as the relative displacements in the transverse direction, the localization 
along the columns, or the vanishing of the tilt modulus are independent of the precise value 
of Uc- However, quantities for which the asymptotic r-dependence is important are crucially 
dependent on Uc- 

The standard determination of Uc in the GVM in the one step case is by minimizing the 
free energy Rr^v in (0) with respect to Uc- This can be done similarly to the case of point 
like disordeilH3. The result for the cosine model in d < 2 is: 
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/3«c = 4T^j.Sr^)/^ (50) 

where ja = (0r(g^ + l)"^- Thus (3uc = 2TE^^^^c^/'^ for d = 1 (ji = 1/2). A similar 
result for a general short range potential can be obtained and involves V'{B). For the cosine 
model if one uses the value of Uc obtained from minimization one gets: 



[UJr, 



^T[^^'/'-{u;1 + J:, + IM)-'/'] (51) 



Expanding the second term in powers of I{uJn) one finds that, unlike what happens in d = 2, 
the linear order does not vanish. As a consequence the solution of (|51|) behaves as /(cUn) ~ u"^ 
at small Un- This solution would lead to an exponential decay of B{t) with r towards its 
asymptotic value. For the classical problem this type of behavior cannot be excluded, though 
it would lead to a surprising difference between d = 2 and d = 1. However, for the associated 
quantum problem there is strong evidence that this solution gives incorrectly the physics, 
for example it leads to unphysical result for the conductivity as we discuss in Section 0. We 
will thus not explore this solution further. 

A look at ( ^41) shows that the general condition for the vanishing of the linear term in 
the expansion in I{ujn), leading to I{ujn) ~ \u!n\ at small Un reads: 

One recognize this condition as the condition of marginal stability, i.e the vanishing of the 
replicon eigenvalue. As discussed in the next section, such a criterion is true at all temper- 
atures. This condition involves the second derivative of the potential and is thus different 
from the equation obtained by minimization of the free energy over Uc- The marginal stabil- 
ity condition is known to arise naturally in the problem of the Langevin dynamics of models 
which have a one step RSB solution in the statics (see discussion in Section [V^) . This is the 
condition that we choose here, on physical grounds. For d = 2 this condition coincides with 
the one obtained by minimizing the free energy studied in section |1V A . 
With this choice of Uc, (3uc = TS^^'^^c^/^ (|5lD becomes: 

I{uJn) = 2Sf [Sr'/' - {ul + Si + IM)-'/'] (53) 



The solution of (|53| ) can be written as 

= T.,f{^c\uJn\/^l) (54) 

where 

2 4 

fix) ~ —;=x when a; ~ (55) 

f{x) = 2 — when x —>■ oo (56) 



Using the above solution and formulae (|41| ) one can compute the correlation functions. For 
fixed r, the a;„ = mode dominates at large x and B grows as: 
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where we have used the value = T / ^/cSi and ^ is the pinning length. 

In the limit of small disorder W , and away from the transition, Ei and thus ^ can be 
computed as a function of the temperature T (for j3 — >• oo). Introducing an ultraviolet cutoff 
Qmax = ^) and rescahng the variables using (|5^) one gets from (|3T|) : 

/■A/VsT daduj 1 T / — 

5 = 2r / ^f,^ . n h ~ - iog(v^^/\/^) (58) 

J [2tt)^ cq^ + cu'^ + 1 + j{y/c\uj\) vrc ^ 



where one has used (^5]). Inserting the value of B ([58|) in ( p8D and using /3mc = T/VcSi one 
gets: 



4T/71 



^ ^ (4Vf//y^)3^2T7^A"^^^^ (59) 

This formula for the weak disorder behavior is similar to the one that was derived for the 
pinning of charge density waves by Fukuyama and LeeS (for T = 0) and extended to take 
into account the quantum fluctuations (here the thermal fluctuations T > 0) by Suzumura 
and FukuyamaS. The formula (^S|) is not identical however to the one of Suzumura and 
Fukuyama since it involves an integration over /(a;) which it determined self-consistently 



in our theory. In Ref. |3g the corresponding I{uj) is zero, which if taken seriously would 
lead to unphysical results for quantities such as C44 or a{uj). In addition, contrarily to their 
solution, it does not contain any unknown prefactor. The scaling with disorder of (|59D 
coincides with the result from the renormalization group analysisil. The phase transition at 
T = Tbg appears on these formula. Additional calculations at finite temperature have been 
performed in Appendix ^ 

At X = 0, -8(0, r) goes to a constant when r 00. This describes the localization of the 
lines due to the correlated disorder. One finds the expression of the localization radius: 

/l = 4-ln(e/«) (60) 

TC^ TTC 

A comparison with the result (^) shows that B{t) grows until r ~ ^ when it saturates at 
Note the difference between the two lengths l± and ^. At the transition one expects that 
^ diverges as ^ ~ exp{b/{TBG — T)") with a = 1/2 from the RG0. Thus relation (60) would 
predict that: 

~ iTsa-T)-r- C'l) 

which could be measured in numerical simulations. 

For the single cosine model in d = 1 and for j3 finite, one does not expect a real phase 
transition (at T > 0) since the system is quasi-one dimensional. The GVM gives however a 
one step RSB solution at finite large enough /3, /3 ~ T. This solution is of the kind discussed 
in Appendix 0: it is discontinuous and disappears when disorder is decreased below a finite 
threshold. It is an artifact of the GVM and should not be interpreted as a genuine phase 
transition in d = 1. However, in the limit (3 = 00 the GVM gives a continuous solution 
which we believe corresponds to a genuine phase transition. 
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C. solution for the general case 



We now present the results for the general case of an elastic manifold in a correlated 
potential. Upon rescaling with (3, the solution for the zero mode a;„ = is identical to the 
problem of uncorrelated disorder in d dimension. However, the values of the breakpoint Uc, 
and the quantities associated to it, [cr](uc) = Si and B{uc) are different. Indeed Si and 
B{uc) are determined by equations (|25|) : as one sees, the difference comes from the second 
equation and B{uc) is determined by fluctuations in d + 1 space rather than d. At T — >■ 
there is no difference and the length ^ is equal to the Larkin-Ovchinnikov pinning length Rf. 
for the problem with point disorder in d dimensionil. However, Si is less renormalized by 
thermal fluctuations and is larger than Si^^^^^orr at finite temperature and thus ^ < -Rc at 
finite temperature. 

Once Si is determined, the equations which determine the connected correlation function. 



(^,0) are closed. The general expression for the correlation function is then: 



B{x, t) = r) - 0(0, 0)]2) = Bd+i{x, r) + Bd{x) (62) 

where the r independent part B^i^x) has formally the same expression as for the li-dim 
problem with point disorder, with different values of Si and Uc- It is the sum of two pieces: 

~ 2T f d'^q ( Si f^dv [a]{u,v) \ 

Bd{x) = -J- / TT—T^ 2^ 2 I V A + / — 2( 2 I r u ^ ~ cos gx 63 

pUcJ (zvr)'* ycg^(cg^ + Si) Jo cq^[cq^ + [a\[UcV)) ) 

This has a well defined limit when /5 — > oo since (5uc goes to a constant. We have dropped 
the term ^ / (^|d^(l — cos(gx)) which corresponds to the contribution of the zero mode 

in the absence of disorder, since it vanishes as /? — > oo. The part -Brf+i(x, r) is truly dW 
dimensional and contains the contribution of all the modes. It reads: 

2T t d^q 1 

This contribution has a different scaling than the previous one. The characteristic length 
scale is ^ = a/c/Si in both x and r, and x ~ r^/^ (diffusion in the localization tube). Bd{x) 



has the same long distance behavior as in the d-dim problem and , if the problem has several 
length scales, exhibits crossover regimes as discussed in Ref. The quantities I{uj), B{uc) 
and Si are determined by equations (^) and (125|). 



Expanding (pOD to lowest order in T gives: 



/(..„) = 2Si[l - (1 + '^n + lM Y^-2y2^ (65) 

^1 

after simplifications using the replicon condition ( P3| ) which gives Si = {—4:j dV" {0)y^^^~'^\ 
(|65|) is solved for any d as: 



/(.;„) = Si/(v/^|a;„|/.^Si) (66) 
where f{z) is determined by inverting the equation: 
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z 



^ = (1 - ^^^i{z)fl^'~'^ - f{z) - 1 (67) 
For d > 2 (|67|) is valid for Un ^ A. The small z expansion of f{z) is: 

Thus the feature that /(a;) ~ |a;| is always verified (with the choice of the replicon condition) 
at T = 0. One can see easily that the above low-T solution will give the correct picture of 
the localized phase. This is because B{t) — B becomes small at large r. To obtain the small 
uj behavior one can expand (PP]): 

I{uJn) = V"{B)^ f dT{l - COs{uJnT)){B{T) ~ B) + A{uJn) (69) 

I Jo 

where A(cu„) = |/o^rfr(l - cos{uj.^t)){H{B{t)) - H{B)) and H{x) = V{x) - xV'{x). The 
part A{uJn) ~ Zu"^ with Z T because the integral: 

Z = -l dT-{H{B{r))-H{B)) (70) 
is convergent at large r as we check self-consistently. Thus one has: 

IM = 2Si[l - (1 + '^n + lM y,-2y2^ ^ ^^^^^ (^^^ 

^1 

where the exact replicon condition Si = {—4:jciV"{B)Y^^^~'^^ has been used. Expanding in 
powers of /(w„) and using A{uJn) ~ Zu"^ at small a;„, one finds: 



, /4(c + Z)Si, , 
^(^«) ~ y 4_j Kl (72) 

This behavior /(cj) ~ Ic^l at small u comes from the decay B — B{t) ~ at large r. 
This property holds for any temperature inside the localized phase, as we confirmed using 
numerical solution of the saddle point equations. 



V. QUANTUM PROBLEMS 
A. models 

The Hamiltonian ([^) introduced in section |I| is also of great interest for the study of 
quantum problems. The dimension denoted as r is viewed as the time and the Hamil- 
tonian becomes the action of a quantum problem. Correlated disorder along the "time" 
direction arises naturally in quantum problems since for these problems disorder is usually 
time-independent. The quantum statistical mechanics of bosons systems corresponds to 
considering r as the imaginary (Matsubara) time and taking periodic boundary conditions 
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along the time direction. The size of the system along the imaginary time direction r is, in 
this representation, the inverse temperature f3 of the quantum problem. The temperature 
T of the corresponding classical problem is a measure of the quantum effects, since it cor- 
responds to h. One of the difficulties in treating quantum problems from the action given 
in section y comes from the fact that all quantities are computed in imaginary time. The 
physical dynamical quantities (such as the conductivity) which are of interest for the quan- 
tum problem come from retarded correlation functions, obtained from the imaginary time 
ones via an analytic continuation. Such a step is in general highly non trivial. Fortunately, 
as we show here, it can be performed directly on the variational formulas, allowing to obtain 
the physical quantities of interest. 

The action (0) corresponds to the general quantum problem of a manifold of internal di- 
mension d {d = corresponds to a quantum particle) in a iV-dimensional random potential. 
One experimental realization is the quantum behavior of flux lines in high-Tc superconduc- 
tors {d = 1, N = 2 for a single line, d = 3, N = 2 for an elastic vortex lattice), which is 
important in low temperature experimentslll. The solution given here for the problem (|I|) is 
exact in the limit N —>■ oo. 

Many quantum problems can also be described by the action (0) and @. It coincides 
with the phase Hamiltonian used to describe charge density wavesil'Q and Wigner crystalJli. 
The density of the system depends on a phase (p by 

p{x) = po cos(Qx + 0(x)) (73) 

where Q is the modulation vector of the charge density wave. The Hamiltonian describing 
the energy cost in variation of the phase is 

where 11 is the conjugate momentum to the phase and M is the mass of the charge density 
wave (usually a large number for phonon- induced charge density waves). If one introduces 
a random potential coupling to density and goes to the action representation one obtains 
exactly the action (P). The current is given by a continuity equation as j ~ dt4>. 

In one dimension (P) also describes the problem of interacting fermions in the presence 
of disorder. The solution presented here describes the Anderson localization of interacting 
particles. Let us consider for simplicity spinless fermions in ci = 1 interacting via the 
Hamiltonian&@ 

n = -tY,clc, + Vj2n,n,^, (75) 

where denotes nearest neighbors. As is well known in one dimension the fermion 

operators can be represented in terms of boson ones0@Ei. The bosons represents charge 
fluctuations of the fermion system. The complete Hamiltonian (^) become^^E^l 



H = — I dx 

2tx 



2 , 1^ 



>ir)(7rH)^ + (-)(a 



(76) 
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n and are canonically conjugate variables and vrll = d^O. The constants v (velocity of 
excitations) and K incorporate all the effects of the interaction V . More generally, the Hamil- 
tonian ( [76D describes correctly all the low energy properties of any one dimensional spinless 
fermion system provided one uses the correct v and K parameters. For the Hamiltonian (|75D 
V and K can be computed from the Bethe-Ansatz solution@ii, but we will consider them 
in the following as parameters. K = 1 corresponds to non-interacting fermions, whereas 
K < 1 (resp. > 1) corresponds to repulsive (resp. attractive) interaction V. Using ( |75D 
one gets for the Lagrangian 

^=7^ [ rfxrfr[-(a.0)2 + vid^(l)f] (77) 
ZnK J V 

In the following calculations we will set f = 1 which can be done by rescaled time and space. 
Therefore in order to identify (^) with (P) one must have: 

c _ 1 
2T ~ 27tK 

The limit K 0, which is the classical limit of the quantum problem (at zero temperature 
(3 = oo), corresponds to the zero temperature limit T — > for the associated classical 
problem in d+1 dimension. One can now introduce a random potential 

Amp = E / dTVit,T)n,{T) (79) 

i 

After replicating the fermions and averaging over the random potential one gets 

£ = -D/rfrE5:<(r)n^(r) (80) 

a,b i 

where (a, b) are the replica indices. In the continuum limit the density operator will have 
Fourier components around q = and q = 2kF, and using the boson representation the 
density becomes 

n(x, t) = -dMx, r) + J_(e'^kFX+2i<p{x,r) ^ ^ /g^x 
TT 2Tca 

where a is a cutoff of the order of the lattice spacing. Physically this means that to describe 
the low energy properties of the system one can only consider the Fourier components of the 
disorder around g = and q = 2kF, and consider them as independent random variables. 
Such an approximation will be valid provided the disorder is weak enough compared to the 
Fermi energy, and therefore if the fermion density is not too small (for a fixed value of the 
disorder). Using the boson representation the disorder term ( pO|) becomes 

J^ = -DJ2 f dxdTdr'^id^rix, r))(9,0'(x, r')) + cos(20"(x, r) - 20^(x, r')) (82) 

As can be shown either on the fermion or the boson representation the g ~ part of the 
disorder can be eliminated by a simple redefinition of variables, as in (|]), which does not 
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(78) 



affect the q ~ 2kF part. Since the current is j = -dt(p such a change has no effect on the 
conductivity, and its effect on other correlation functions can be simply computedii. In the 
following we therefore keep only the q ~ 2kF part, which can lead to the localization of 
the fermions. The total action we consider in the following is identical to model (P) with 
W/i2T) = 2D/{2Tiafv. 

A similar mapping can be made for repulsive bosons in one space dimensionSEa. The 
resulting Lagrangian is identical to ( |77D with the identifications v/{7iK) = po/m, and v"^ = 
l/^npom). K is the compressibility, po the average density, and m the mass of the bosons. 
The excited states in the absence of disorder are sound waves with phase velocity v , which 
are the phonon modes typical of a Bose superfluid. The existence of such modes is sufficient 
for true superfluidity to existS. 

The point K = 1 is the best known. It corresponds to non interacting fermions, or 
equivalently to bosons with a hard core repulsion. There, it is well known that all states are 
localized by disorder, and that the conductivity behavesEj-ea as (j{oj) ~ a;^log^(a;). For small 
disorder one can use the renormalization group (RG) and one finds a transition between a 
regime where the disorder is irrelevant, and a regime where the disorder renormalizes to large 
values (there is no perturbative fixed point accessible in the localized phase). The localized 
regime is found to happen for K < 3/2, i.e T < Tbg = 37rc/2 with the identification (fTHI). 
This coincides with the value predicted by the variational method (see section IV B| and 
Appendix A). The localization length, the crossover towards the localized regime, and the 
high-frequency part of the conductivity can be also computed using RCi. When ^ the 
hamiltonian becomes similar to the one describing the pinning of a charge density wave, for 
which a self-consistent formula for the pinning length was obtained by Fukuyama and Lee. 
The conductivity was computed analytically for K = using a transfer matrix approach 
and is also found^^ to behave as (y{uo) ~ uo"^ \og^ {uo) ^ as for K = 1. 

Most of the physical properties can be derived from the solution of the GVM in imagi- 
nary time obtained in section |1V B\ Contrarily to the case of the elastic manifold, the length 
^, e.g given by (|59|) (the length at which the phase disorders) is interpreted as the local- 
ization length for the interacting fermions. This is because it set the frequency scale in the 
conductivity as can be seen on the GVM result below and on the RG. For K = 1 it coincides 
with the usual definition of the localization length at the Fermi energy. At a given point in 
space formula ( |60D and (^) show that the phase fluctuations saturate to a constant at large 
time. Thus, in the localized phase, at each point there is a well defined average value of the 
phase that can be viewed as constant in time but which adjusts spatially in order to take 
advantage of the random potential. This agrees with the physical picture of Fukuyama and 
Leeii. 



B. conductivity: generalities 

In addition to the results which can be directly transposed from section let us examine 
here dynamical correlation functions. The most interesting quantity is the conductivity. We 
will compute it using the Kubo formula (see Ref. Ch. 3 for a derivation, and also Ref. 
g): 



21 



1 yj^ 

(^{(l:^) = -{ + Xret(g,^)) (83) 

UJ TC 

Here u is a. real frequency associated to real time t. It is understood that u always stands 
for u + i6, i.e contains an infinitesimal imaginary part. The first term is the diamagnetic 
contribution and the second term is the retarded current-current correlation function: 

Xretiq,co) = -il dte'-'9it){[j^iq,t),jiq,0)]) (84) 



The conductivity is defined as o"(c(j) = limq^oa{uj,q). For the actual calculation one must 
work with imaginary time r = it, use the imaginary Matsubara frequencies iUn and then 
perform the analytic continuation iun ^ uj + i6. Introducing: 



X{q,iu;n) = - rrfre*""^(T.J+(g,r)j(g,0)) (85) 



et done XreM,^) = ^^n)|ja.„^c^+j5- Using j = ^dtcj), the integration by part over time 
in (^) cancels the diamagnetic termcHl and one finally finds the general formula for the 
conductivity: 



^(^) = HiTMq = 0, u;„)0(g = 0, -cu„)) , (86) 

One can apply this formula to a system of interacting spinless fermions (or bosons) in the 
absence of disorder. The correlator is then given by (Tr0(g, a;„)0(— g, u;„)) = T/{cu!^ + cg^) 
and thus one finds: 

a{uj) = hm „, ' (—^^^] = = D (6{uj) + -P{-)) (87) 

where P denotes the principal part. D = T/iTcc) = vK is the charge stiffness and the 
strength of the Drude peak. The charge stiffness D can be defined quite generallylll by 

(Treg(cu). When D > it means that the system is a perfect conductor. 
To explicit the content of the flux line-quantum boson analogy it is useful to indicate 
the relation between the conductivity of the quantum particles and the tilt modulus of the 
corresponding classical system. Using (EBI) and (BB) one has: 



T 



C44:{uJn) 



-TT^iujcx{uj) (88) 



This relation between C44, and a dynamic quantity, (7(0;), involves a non trivial analytic 
continuation and one should be careful before drawing conclusions from one system to the 
other. Taking the limit — > in the above equation one gets: 

T 

— = lim iT^Im(uJ(r(uj)) = nD (89) 
C44 '^-*o 

where the last relation follows from the Kramers-Kronig relations. Thus the tilt modulus C44 
is related to the charge stiffness D of the quantum system. Both D and C44 correspond to 
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"dynamical" quantities, i.e the hmit g — is performed before the c<j — > 0. Note that both 
quantities are unrenormahzed by interactions for a gahlean invariant system. In general, 
I/C44 ~ D should be distinguished from "static" quantities such as the superfluid density 
Ps ~ d'^F/ which can also be obtained from the current-current correlation function, but 
in the opposite limit 00 first then g — ^ (F is the free energy and $ an external flux). A 
priori D and ps are different quantitieJH (for free fermions Ps = and D 0). At /3 = 00, 
D ~ d'^Eo/d^^ where Eq is the ground state energyiil. In = 2 and P = 00 a relation 
between C44 and ps applied^ but only when cqq = and thus only in the liquid phaseiH, i.e 
in the superfluid phase for the bosons. 

In the presence of disorder, in the Bose glass phase, both o"(ti; = 0) and C44 are expected 
to vanish. A point which deserves further investigation is whether the two quantities D and 
Ps may generally become proportional in the presence of disorder. 



C. conductivity computed from the GVM 

In order to compute the conductivity in the most general case one needs the analytic 
continuation of ( p9| - pl| ). Let us first examine the limit — 0, i.e T — 0, where the 
continuation is straightforward. The analytical continuation to real frequency u can be 
done directly on the equation p^ ) and gives: 

H^) - -4r'(0) / Jg, - ,^.^^,_^,,.^,(,) ) (90) 

In the case of the single cosine model (i.e for interacting fermions) in d = 1 it gives: 

liu) = 2Ef 'p^'/' - (Si -co^ + Huj))-'/^] (91) 
The solution of this equation, and of the general case of an arbitrary V{B), can be written: 

= Si/(v^za7/v^) (92) 



where the scaling function f{z) was computed in Section IV. 
The expression (^) for the conductivity can be rewritten as: 

-Tu; I"{u;)+z{I'{u;)-cu;') 
vr^ /"(^)2 + (/'(^)-c..2)2 ^y^J 

where I{iuJn) — ^ lioj + i5) = + il"{uj). It can also be written in the scaling form: 

^(^) = ^siV^zu;/^,) (94) 

where s{z) is an analytic scaling function of z. We are primarily interested in the small 
frequency behavior. Writing f{z) = 1 + bz + O^z"^), with b = 2/^/3 for d = 1, one finds the 
asymptotic behavior at small frequency: 
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Re{a{uj)) 
Im{a{ijj)) 



bTu^ 

7r2Si V Si 
-Tiu 



(95) 



To compute the conductivity at positive temperature T > is more difficult because the 
analytic continuation of (29-^) becomes non trivial. The continuation of the mass term 
Si(l — Snfl) give simply Si, but the continuation of the self-consistent equation for I{iUn) 
is more involved and is performed in appendix 0. The result reads, for the single cosine 
model: 



rret / 



AWe 



-2B 



T 



/3 



dr i e 



= ^ /J"(iMyl(«){e-"l^l+2iV„cosh(ur)) 



+2 



dte'^^'lm 



^ /^°°d«A(u)(e-^"«+2Af„cos{ui)) 



(96) 



where iV„ 



l/(e'^" — 1) is the Bose factor and A{u) is the q integrated spectral function 
d'^q I"{u) 



A{u) 



{2tiY (g2 - m2 + Si + r{u)f + {r{u)y 



(97) 



with the property A{—u) = —A{u). One must check that indeed in the limit of -f^ 
(equivalently T 0), when one can simply expand the exponential one recovers the above 
expression (|90D. The expansion gives: 



JM = IQWe 



duAiu)-{- + 



(98) 



Using again ( p5| ) one obtains back (pO|). 

We do not expect that finite T {K > 0) will change qualitatively the low-frequency result 
apart from a change in the coefficient b when T varies, presumably an increase of b with 
T. The high frequency behaviour however is strongly dependent on T and is computed in 
Appendix O. 



D. discussion 

We have found a conductivity cr{uj) ~ u'^ in any dimension d using the GVM. This is 
a direct consequence of the iu term being generated in the self-energy. While this term 
appears naturally when there is full RSB, in d = 1 for the single cosine model, it appears 
only as the consequence of the marginality condition in the replica language. The other 
choice (minimization over Uc) leads to I{uj) ~ and to a gap. Indeed generically, simply 
adding a mass term to a uj"^ term in the self-energy leads to a gap in the excitation spectrum. 
If one has {Tr(p{q,u!n)(p{—q,uJn)) = T/{cq^ + du"^ + Si), then using (|86| ) one finds: 

^(^) ~ V "^'2 -A ~ - \/^) + ^n^^) (99) 
1^1 — uj-^ — 10 ^ a; — Li 
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Thus there is a gap and no absorption at small frequency. 

Within the framework of the variational method the present result, ~ tu^, is ob- 

tained by taking into account only small oscillations around the equilibrium position (in the 
following we call "phonons" such excitations, not to be confused with real phonons which 
are absent here). Our result is in reasonable agreement with previous results in = 1, up 
to logarithmic corrections. Previous results were: (i) for free fermions at i^' = 1, cluj) ~ u"^ 
by Mott, and cr{uj) ~ cj^log^a; through exact calculation by Berezinskii@@ and through a 
physical argument by Halperin, (ii) in the classical limit K = 0, a{uj) ~ cj^log^u; using a 
transfer matrix techniqueEj by Feigelman and Vinokur (FV). In Mott's physical derivation, 
the conductivity is estimated from computing the energy absorbed by the system of elec- 
trons. An external driving frequency u will cause electrons within u of the Fermi surface 
to make transitions from one localized state to another. The density of the active electrons 
is thus ~ u and since they absorb a photon of energy u, one finally finds a conductivity 
~ uj"^. Thus in this derivation, it seems that the important excitations are, when cast 
back in the phase language, the soliton-like excitations of the phase, corresponding to an 
electron being transferred from one localized state to another. In FV derivation however, 
only "phonons" modes are considered. 

One may thus speculate that both solitons and "phonons", i.e small oscillations of the 
phase, give the same contribution u;^log^(u;) to the conductivity. One must note that since 
we have considered Gaussian disorder, corresponding to the limit of a large number of very 
weak pins, it is always possible to create phonon-like excitations down to = 0. For 
Poissonian disorder, on the other hand, phonon-like excitations should be impossible^ for 
frequencies lower than the strength of an individual pin V. In that case the d = 1 wire 
is cut into segments with fixed boundary conditions, and the conductivity from "phonons" 
is exponentially activatedBil. On the other hand soliton like excitations are still allowed. 
One can therefore expect that above uj = V both phonons and solitons will contribute to 
the ci;^log^(a;) term in the conductivity, while below u = V, only the contribution of the 
soliton part should remain. In the near classical limit, for small K, solitons are energetically 
much less favorable than phonons and the contribution of the soliton part should therefore 
be small. The conductivity should therefore remain as iu'^\og^{iu) down to arbitrary low 
temperatures but the prefactor dropping dramatically below uj = V. In the other case of 
large quantum fluctuations, for example free electrons, there should be practically no effect 
of Poissonian disorder. 



VI. CONCLUSIONS 

In this paper we have extended the application of the gaussian variational method (GVM) 
to quantum systems with disorder and, using the analogy between the two problems, to 
classical elastic systems with correlated disorder. By construction this method is exact in 
the limit of a large number of components, i.e in large embedding dimension. It allows to 
describe the equilibrium statistical mechanics of both problems. A glass phase with replica 
symmetry breaking (RSB) is found at low temperature. 

Compared to the case of point disorder treated in Refs. the solution in the case of 



correlated disorder presents three new non-trivial features. These are not merely technical 
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complications but originate from the physics being very different in the present case. 

First a correct treatment of the boundary conditions is crucial. The disorder being 
correlated along r (imaginary time for the quantum particles) a single mode, the mode u = 0, 
controls most (but not all) of the physics in the glass phase. The RSB is confined to this 
mode and there is a non trivial mechanism to generate an effective mass term for all the other 
modes, leading to localization by disorder. The associated physics which naturally emerges 
from the GVM is that of small distortions around a classical background configuration. 
The method therefore confirms the physical picture developedH by Fukuyama and Lee. In 
addition it allows for a calculation from first principles of all correlation functions. The 
properties of the classical configuration is fully described by the RSB and is not assumed 
from the start. For classical elastic manifolds we have obtained the pinning length ^, the 
localization length l± and the non trivial frequency dependence of the self-energy I{uj). 
The GVM predicts a space vs "time" scaling x ~ with z = 2 within the localized 
phase and up to the transition. For the single cosine model, this agrees with a recent RG 
predictionii in d = 4 — e obtained directly from the RG fixed point for the mode u = 



previously derived in Ref. |2^,|29|. These results could in principle be checked experimentally 
by measuring correlations along the columns using neutron diffraction. We have computed 
the tilt modulus C44 and shown that C44 = cxd in the Rose Glass phase, a signature of the 
transverse Meissner effect predicted in Ref. ^. By including higher harmonics and realistic 
elastic energy the GVM can be used to obtain a more detailed description of the Bose Glass 
phase. 

In order to apply the GVM to quantum problem and compute dynamical linear response 
quantities, such as the conductivity, we had to perform a non trivial analytic continuation to 
real time, t = it. We showed that within the GVM such a continuation can be performed in 
full generality. We computed the frequency-dependent conductivity (j{uj) for phase Hamil- 
tonians, relevant e.g. for Wigner crystals and Charge Density Waves. These systems are 
localized by the disorder and we computed the localization length ^. We find that o"(u;) ~ tu^ 
at low frequency in any dimension. We have also studied interacting fermions and bosons 
m d = 1. The method then predicts the localization-delocalization transition correctly. It 
yields an analytic form for cr(co') which is in agreement with known high frequency results 
and which at low frequency also agrees, up to logarithmic factors, with the exactly solvable 
cases (i.e non-interacting fermions and classical limit). Such a behavior of the conductivity 
results from a new "friction" term iuj generated by the random potential as predicted by 
the GVM. Finally, a general relation between the tilt modulus of the classical problem and 
the conductivity of the quantum problem has been derived. 

The third new feature arises when looking at the d = 1 quantum problem, modeled by a 
single cosine Hamiltonian. We found that the natural way to obtain a physically correct dy- 
namical response (conductivity) was to use the so-called marginality condition at the replica 



saddle point. It is known that the marginality condition arises naturally in the solutionE^I of 
off-equilibrium Langevin Dynamics (LD) for spin glass models. The equilibrium quantum 
statistical mechanics studied here thus plays the same role with respect to the full problem 
of quantum dynamics as the equilibrium FDT-regime does in the off-equilibrium Langevin 
Dynamics (LD) of a classical system (in which the fluctuation dissipation theorem (FDT) 
and time translational invariance holds, and thus called the FDT-regime). This scenario 
agrees with the fact that the quantum equilibrium dynamical equations that we obtained 
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( p9| - |3T|) depend on the quantity Si only. This quantity is determined at the breakpoint 
from the RS side (short scales) independently of the details of the RSB at larger scales. 
This is completely analogous to the role played by the "anomaly" M with respect to the 
FDT-regime of the off-equilibrium LD. The LD for the d-dimensional classical problem (i.e 
the = mode of the quantum problem) was studied in Refs. ^^"^ expressions 

for the anomaly M and for the large time limit hi of the displacement in the FDT-regime 
were derived. One can check explicitly that in the limit T 0, Si and B in the quantum 



problem are identical to the anomaly and hi from Refs. |58| , |57| . Thus we speculate that there 
also exists an off-equilibrium quantum aging regime, at finite waiting time tw The regime 
studied here of equilibrium time-translational invariant quantum dynamics being reached 
only at infinite waiting time — oo and fixed t = it. This aging regime might be relevant 
for computing some response functions, in particular the nonlinear responses. 

As we have discussed, the GVM studied here takes into account only the small gaussian 
deviations around a non trivial, spatially dependent, equilibrium (i.e static) configuration. It 
would very interesting to extend the method and take also include kinks and solitons. This 
would allow to compute the response to a finite transverse field (i.e study the transverse 
Meissner phase) as well as non linear dynamical quantum responses and variable range 
hopping phenomenon. 
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APPENDIX A: REPLICA SYMMETRIC SOLUTION 

Assuming replica symmetry in (|TBp one finds: 

G-^ = c{q^ + u:l)+'^ f dr{l-cos{u^r))[V\B{r))- f duV\B))] 
I Jo Jo 
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a{uJn) = -^6n,oV'{B) (Al) 



T 



where 



2T r d'^q 
P (2vr)' 

B{r) = y EX - ^osM) (A2) 

Since the a; = mode is infrared divergent in q, B = oo for d < 2. The equations simplify 
in: 

= c{q' + ujI) + ^ dT{l - cos(.;„r))V^'(S(r)) 
2T d'^a 

B{r) = y EX - cos(^nr)) (A3) 
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Contrarily to the RSB case no good T = limit exists. Thus the solution cannot describe 
a low temperature localized phase with a finite localization length. A localized solution of 
( |A3| ) can be found numerically (adding a mass term), but has no good limit /3 — > oo since 
the integral in the first equation produces a term proportional to (3. 

The RS solution describes however the high temperature phase. For power-law models 
such a phase exists when in the short-range or marginal case, 7 > 2/(2 — d). For the single 
cosine model it exists in d = 1. It also exists in ci = 2 but since Tc cc (3 it disappears in the 
limit /9 — *■ cxo. Starting from the elastic theory of the fiux lattice in ci = 2 correlated disorder 
always lead to the Bose glass phase. Only the liquid, such as the case c? = 1 of lines confined 
to a plane, may show a transition. 

In the high temperature phase, the tilt modulus C44 is renormalized, and given by (in the 
limit (3 = 00): 

1 (It 

C44 = c+-y^ -rV(i?(r)) (A4) 

The divergence of the integral to first order in perturbation theory signals the transition at 
Tbg- Foi' the model which describes the vortex lattice one finds V{x) = —Wexp{—KQx/2), 
Tbg = Qttc/Kq. Within the self consistent method one has: 

cu = c+—^ r2exp(--=^log(r/A)) (A5) 

I Jo ZTT v47rC44C 

there is no divergence of C44 at the transition, but the solution for C44 ceases to exist below 
Tbg- One recovers correctly some of features of the KG approach. 



APPENDIX B: STUDY OF THE RSB SOLUTION FOR THE SINGLE COSINE 

MODEL IN D = 1 AT FINITE T 

We now study the equations for the single cosine model in ci = 1 at finite T. We will 
study the limit where Si/A <C 1, i.e ^ S> a. This happens in two cases: (i) at weak disorder 

— > at fixed T (ii) near the transition T —>■ Tbg- We will only study the case /3 = 00. 
We take c = 1. 

The closed system of equation (for the displacements, self energy and effective mass) is: 



B(r)=T 



^ duj cos(ujt) 



I{uj) = / cir(l - cos(^r))(e2^M - 1) (Bl) 



T 

(Si)3/2 = 4iye-2^(°) 

where = A is the ultraviolet cutoff, and we have denoted by -B(r) = B{t) — B. All 
variables can be rescaled when Si — >• 0. We thus define r = (Ei)~-^/^ti, oj = (Si)^/^ty, and 
as before I (to) = Si/(w). The cutoff becomes Wm = (Si)~^/^f2 and we are thus interested 
in the limit Wm — ^ 00. Eqs. (pT|) become: 
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B{u) = T j_ 



27r 



oo 

TJo 



f{w) = ^ 1' du{l - cos(wM))(e2^(") - 1) 



= 4l^e-2^(°) w;^ = {^^)-'/^Q (B2) 

The first two equations have a well defined limit when Wm ^ oo as long as T < Tbg = 37r/2. 
Thus the functions B{u) for u > and f{w) for w < oo also have a well defined limit. The 
approach of the BG transition from below is characterized by divergences in the integrals 
at small rescaled time u. As we show below, when T < Tbg one has f{w) <C w'^ at large 
frequency ^ 1. Thus one has B{u) ~ —{T/tc) \ia{u/u*) at small u, and reporting in the 
second equation one gets: 

f{w)r^— I dui — )~ il — cosiwu)) (B3) 
T JO u 

One sees that this integral converges at small u provided T < Tbg = 37r/2. There is thus a 
well defined weak disorder continuum limit in c? = 1. 

Analyzing further the solution one finds the following. For T < 7r/2 (i.e K < 1/2 ) 
lim^^oo fi.'uj) = /(oo) is finite. For 7i/2 < T < 37r/2 one has f{w) ~ w°' with a = — — 1 
(for 1 <^ w <^ Wm- One sees that the point K = 1 (free fermions) corresponds to f{w) ~ w 
both at small and large frequency. The transition at T = Tbg = 37r/2 occurs when a = 2, 
i.e when the term /(a;) becomes of the same order as cu^ at large frequencies. Of course the 
small frequency behavior is always the same ~ |w| as discussed in the text. 

Using the formula (|93|) for the conductivity and the above result I{uj) ~ uj"^^'^ at large 
uj, one finds: 

This power-law behavior of the conductivity coincides with previously known resultsiiS. 
The GVM therefore describes the correct physical behavior for the conductivity for all 
frequency scales. 



APPENDIX C: REMARKS ON THE VARIATIONAL METHOD 

For finite number of components N the Hartree or Gaussian Variational Method is only 
an approximation. It is important to determine whether it is a good approximation. It is 
in principle possible to answer in a expansion, but this is a rather difficult task. The 
GVM seems to work well e.g for the cosine model in the case of point disorder where it can 
be compared with RG calculations. In d = 4 — e the values of the amplitude of logarithmic 
growth of displacements was found^ to agree within 10% in the FRG and GVM. In = 2 
it predicts the transition temperature Tq exactly. 

In this Appendix we illustrate some of the artifacts produced by the method. When 
properly understood these examples may help decide to which extent and in which cases the 
results of the GVM can be trusted. The following simple model is studied: 
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j d'^x^{V(f)y - g cos (l){x) (CI) 

The variational Hamiltonian is i^o = / d'^q/ (27r)'^)(l/2)G~^(g)0(g)0(— q') and the variational 
free energy F^ar — —T\nTrex.p{—Ho/T) + {H — Hq)h^ is found to be: 

F,ar T f d'^q ^^^Tc f d^q , T f d^q 

IT = - 2 y (2^ ^"^^^^^^^ + T y (2^^ ^^^^ - ^exp(-- J j^,G{q)) (C2) 

where V is the volume of the system. The saddle point equation thus gives G~^{q) = cq'^ + a 
with: 

T f d^q 1 , 
a = ,exp(--y^-^) (C3) 

Let us consider = 1, i.e a string directed along z, in a periodic potential independent 
of z. In the limit of a large ultraviolet cutoff, qmax = A — > oo, one finds, up to a diverging 
but cr-independent contribution: 

^'-(x-^e-V-) (C4) 



y 16c 



w 



2\ 



rhcrc we have defined x = i^/ca/T and g — 16cg/T^. The total field fluctuation is (0 
l/2a;. One sees that at small g, F^ar is a monotonic function with a minimum at a; = 0. 
The first term, the entropy, dominates and the string is delocalized, (0^) = oo, which is the 
correct result. However, for g > e^/4 a second minimum at finite x appears (it appears at 
X = 1/2 for ^ = e^/4). There are then three roots to the saddle point equation. The absolute 
minimum remains in a; = until g reaches g — e (when x — 1), at which the (unphysical) 
finite X solution becomes of lower free energy. 

The fact that in d = 1 the GVM predicts a transition at ^ = e to a localized phase is 
an artifact. The physical reason for this result is very clear, however. The first term in the 
variational free energy is the entropy cost ~ T^/ci?^ of confining the string to a transverse 
extension R, i.e the kinetic energy of a bound state in the quantum mechanical analogy. 
This is a correct estimate. The energy cost of delocalizing the string away from a minimum 
is, however, strongly overestimated at large g. Indeed it is computed by assuming that the 
localized packet is gaussian, while in fact the packet will be strongly peaked near the minima 
of the potential. The true delocalization will thus proceed via kinks and be less costly in 
energy. At small g this effect is not a problem, and one expects the method to work well. 

In the case d = a similar artifact occurs, though in a weaker sense. The equation 
a = g exp{—T / {2a)) has always the physically correct solution a = 0. If g/T > e/2 two 
additional solutions appear: a — a+ and a — a-, such that at g/T — e/2, a- = a+ = T/2. 
However in the GVM in d = entropy always wins and the absolute minimum in the 
variational free energy is at (0^) = oo, which is the correct physical result. At large enough 
g, however, a second unphysical minimum appears, which leads to some artifacts in the 
presence of an extra mass term. 

Our conclusion is thus that one may trust the GVM only for weak coupling or disorder. 
A transition at intermediate couplings, in low dimension, which cannot be tuned to zero, is 
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usually the sign of an artifact. We believe a similar artifact to occur in the case of d = 1 
with point disorder (see Appendix C in Ref. 2£ and Section IV B| above). This artifact 
may be indicative of tendencies of the system, e.g. to be more glassy, but not taken as the 
sign of a true transition. In the disordered case, one may also be able to determine when 
statistical entropy and disorder energy are incorrectly evaluated. One generally expects that 
the main corrections should be produced by kinks or instanton configurations, while small 
displacements, e.g. phonons, should be described correctly. 



APPENDIX D: ANALYTIC CONTINUATION 



We illustrate the method of analytic continuation to real time on the single cosine model 
(^), with V{x) = —Wexp{—2x), which is relevant for interacting fermions in ci = 1. The 
extension to arbitrary V{x) is straightforward. 

We want to perform the analytical continuation of: 

F(icj„) = 3ll| j dr(l - cosKr)) ^ '^^^ _ ^ j 

where I{iuJn) is the self energy with Matsubara (imaginary) frequencies. Physical quantities 
are expressed in terms of the retarded self-energy /(cu) = I'{oj) + il"{uj) obtained by the 
continuation iUn uj + i6. The term coming from the 1 in the integral ( pi| ) is a simple 
constant term, independent of frequency and is therefore unchanged under the analytic 
continuation. Difficulties come from the cosine term. This term is the Fourier transform in 
Matsubara frequencies of the imaginary time function 

/cos(r) = e^ ■' c(.^+,^)+Ei+7(,.p) _ ^ P2) 

since the Un = 27171/(3 are Matsubara boson frequencies. In general, obtaining r^^{t) from 
(P^), is a difficult task. Here, however, I is computed using the Gaussian variational ap- 
proximation. Within this approximation Jcos is itself a Green function. One can therefore 
use for /cos the well known relation between the real time time ordered boson Green function 



Icos{t)=9{t)hit) + 9i-t)l2{t) (D3) 

and the retarded one to get 

to perform the analytical continuation, it is therefore sufficient to know the real time, time 
ordered, function Icosit). 

In order to do so, we use a spectral representation. 

1 ~1 r°° 1 
Gc[q,iu;p] = = — / duA{q,u)- (D5) 



where the spectral function A{q,u) is the imaginary part of the retarded function (see Ref. 
50| sections 3-2 to 3-5) defined by 
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A{q, u) = Im{Gc[q, iuip ^ ui + i8\) 



r>) 



(D6) 



The above identity uses 5(^7 — n) = ^/m[l/(u; + i5 — u)]. The real part of (P5| ) gives 
i?e((j'c[Q'; "^^p ~^ uj + i5]) = — J duA{q, u)P{l/ (u — u)), which is the Kramers-Kronig relation. 
The term in the exponential in (p2|) then becomes 



-4T f d'^q 



p7T J {2tT 



u 





iujp — u 



(D7) 



Using the fact that A[q, —u) = —A[q, u) for bosons systems one can rewrite (|D7|) as 
AT f d'^q , ^ , . 2m 



/(2^r^"? 



cosi^c^pTj^- — -A{q, u) 



(371 J (27r) 

One recognize in (^^) the Fourier transform of a simple phonon propagator which gives 
4T f d'^q roo 



(D8) 



TT J (27r)"' Jo 



rfn[e-l"l" + 2N^ cosh(Mr)] u) 



(D9) 



where Nu = l/(e^" — 1) is the boson occupation factor. The time continuation to real time 
r = is now easily performedll and gives 



4T f d'^q 



TT J {2ny Jo 
Using (|D3| ) and (P4|) one finally obtains 



rfM[e-^l*l" + 2Ar„cos(ut)]A(g, 



u 



(DIO) 



AWe 



T 



+2 / dte'^^'lm 



dr (^e^ /7d«A(M)(e-"l-l+27V„cosh(«r)) 
duA{u){e-'^^+2NuCosiut)) 



(DU) 



where A{u) = J ■7^:4aA{q, u) is the q integrated spectral function. 



(2^) 
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